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Exact solutions of the gravitational field equations for a mixture of a null charged strange quark 
fluid and radiation are obtained in a Vaidya space-time. The conditions for the formation of a naked 
singularity are analyzed by considering the behavior of radial geodesies originating from the central 
singularity. 



I. INTRODUCTION 

O; 

It is generally believed that strange quark matter, consisting of u, d and s quarks, is the most energetically favorable 
state of baryon matter. Witten [1] suggested that there are two ways of formation of the strange matter: the quark- 
hadron phase transition in the early universe and the conversion of neutron stars into strange ones at ultrahigh 
densities. Some theories of strong interactions, e.g. quark bag models, suppose that the breaking of physical vacuum 
takes place inside hadrons. As a result the vacuum energy densities inside and outside a hadron become essentially 
i different and the vacuum pressure on a bag wall equilibrates the pressure of quarks thus stabilizing the system. 
£NJ ' There are several proposed mechanisms for the formation of quark stars. Quark stars are expected to form during 
the collapse of the core of a massive star after the supernova explosion as a result of a first or second order phase 
~~ * ■ transition, resulting in deconfined quark matter [2]. The proto-neutron star core or the neutron star core is a favorable 
environment for the conversion of ordinary matter to strange quark matter [3] . Another possibility for strange star 
formation is that some neutron stars in low-mass X-ray binaries can accrete sufficient mass to undergo a phase 
transition to become strange star [4] . Such a possibility is also supported by an unusual hard X-ray burster [5] . This 
i mechanism has been proposed as a source of radiation emission for cosmological 7-ray bursts [6]. Other possible 
astrophysical phenomena related to strange stars are discusssed in Ref. [2] . Hence the problem of the collapse of the 
-H , strange matter is of much interest in the study of strange quark star formation, dynamics and evolution. On the other 

hand the theoretical understanding of the collapse is of fundamental importance in general relativity. 
q I Usually it is considered that quark matter is formed from a Fermi gas of 3^4 quarks constituting a single color- 
singlet baryon with baryon number A. The theory of the equation of state of strange matter is directly based on the 
fundamental QCD Lagrangian [7] 
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. £j i L QCD = l £ F%F^ + (*fd* - vfA'Jj - "»/) 1>- (!) 

The subscript / denotes the various quark flavors u, d, s etc. and the nonlinear gluon field strength is given by [7] 

F% = MS - d » A l + gfabcA^Al. (2) 

QCD predicts a weakening of the quark-quark interaction at short distances (or high momenta Q 2 ), because the 
one- loop series for the gluon propagator yields a running coupling constant. 

Neglecting quark masses in first order perturbation theory, the equation of state for zero temperature quark matter 
is [1,7] 

P =l(p-m, (3) 

where B is the difference between the energy density of the perturbative and non-perturbative QCD vacuum (the 
bag constant) and p,p are the energy density and thermodynamic pressure of the quark matter, respectively. Eq. 
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(|3j) is essentially the equation of state of a gas of massless particles with corrections due to the QCD trace anomaly 
and perturbative interactions. The vacuum pressure B, which holds quark matter together, is a simple model for the 
long-range, confining interactions in QCD. At the surface of the star, as p — > 0, we have p — > 4B. The typical value 
of the bag constant is of the order B — 57MeV/ fm 3 w 10 15 g/cm 3 [2]. After the neutron matter-quark matter phase 
transition (which is supposed to take place in the dense core of neutron stars) the energy density of strange matter is 
p ss 5 x 10 14 g/cm 3 [1]. Therefore quark matter always satisfy the condition p > 0. 

It is the purpose of the present Letter to obtain an exact non-static solution of the Einstein field equations for 
a collapsing charged null strange fluid in a Vaidya space-time and to study some of its singularity properties. The 
Vaidya geometry permitting the incorporations of the effects of radiation offers a more realistic background than static 
geometries, where all back reaction is ignored. The exact solution obtained represents the generalization, for strange 
matter, of the collapsing solutions previously obtained by Vaidya [8], Bonnor and Vaidya [9], Lake and Zannias [10] 
and Husain [11]. For a general method of obtaining spherically symmetric solutions in Vaidya geometry, see Ref. [12]. 
The structure and properties of singularities in the gravitational collapse in Vaidya space-times have been analyzed 
in Refs. [13-17]. 

The present Letter is organized as follows. In Section II the gravitational field equations are written down and the 
general solution is obtained. The conditions of formation of a naked singularity are analyzed in Section III. In Section 
IV we discuss and conclude our results. 



II. SPHERICAL COLLAPSE OF THE STRANGE QUARK NULL FLUID 

In ingoing Bondi coordinates (it, r, 9, tp) and with advanced Eddington time coordinate u = t + r (with r > the 
radial coordinate and r decreasing towards the future) the line element describing the radial collapse of a coherent 
stream of charged strange matter is 



ds 2 



^ 2m (u, r) 



du 2 + 2dudr + r 2 (dO 2 + sin 2 6dip 2 ) . (4) 



m (u, r) is the mass function and gives the gravitational mass within a given radius r. 

Equilibrium configuration of masless strange quark matter has equal numbers of u, d and s quarks and is electrically 
neutral. If the mass of the strange quark is not zero, strange quarks are depleted and the system develops a net positive 
charge. Since stars in their lowest energy state are supposed to be charge neutral, electrons must balance the net 
positive quark charge in strange matter and the condition of charge neutrality requires |n u — ^n^ — -|n s — n e = 0, 
where n is the particle number density [2]. Hence quark matter is usually neutral. In principle, rather extreme 
astrophysical conditions could lead to a charged strange quark matter astrophysical configuration and, for the sake of 
generality, we shall also consider this possibility. 

We shall write the matter energy-momentum tensor in the form [11,12] 

where 

T^>=p:(u,r)y v , (6) 
is the component of the matter field that moves along the null hypersurfaces u — const., 

= {p + p) (l„n„ + l v n„) + pg^, (7) 
represents the energy-momentum tensor of the strange quark matter and 

E nv = TZ \ F m F v - -79^F al3 F af3 ) , (8) 



2m(ti,r) 



4tt V 4- 

is the electromagnetic contribution, and n M are two null-vectors given by = 69$. and = | 

5qI, so that l a l a = n a n a = and l a n a = —1 [11,12] (with 8$ the Kronecker symbol). The energy density and 
pressure in Eq. (Q) have been obtained by diagonalizing the energy-momentum tensor obtained from the metric [11]. 
The electromagnetic tensor F^ v obeys the Maxwell equations 
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(9) 
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Without any loss of generality we take the vector potential to be [10] 
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g(*0 f(M) 



(10) 



(11) 



with g(it) being an arbitrary function. From the Maxwell Eqs. (|9j), (10) it follows that the only non- vanishing 
components of F^ v are F ru = —F ur = and 



^ = ^^(-1,1,-1,1). 



For the energy-momentum tensor (^) the gravitational field equations take the form 

1 dm (it, r) 



r 2 du 



4:77/1 (it, r) , 



(12) 



(13) 



2 dm(u, r) „ , . Q 2 (u) 



dr 



(14) 



ld 2 m(u,r) o , . g 2 (u) 
r or" 1 r 4 



(15) 



By using the bag equation of state of strange matter given by Eq. (^) , we obtain the following equation describing 
the dynamics of a mixture of null fluid, strange matter and electric field: 



d 2 m(u,r) 2 dm(u,r) 2q 2 (u) 32nB 



dr 2 



3r dr 3 r 3 



By means of the substitution 



m (u, r) = mo (u, r) H — r — 



3 2r ' 

we obtain a new unknown function mo (it, r), which satisfies the equation 

d 2 mo(u,r) 2 dm (u,r) 



dr 2 



3r <9r 



having the general solution 



m («,r) =C(m) r 1/3 + F(u), 



(16) 



(17) 



(18) 



(19) 



with C (u) and F(u) being two arbitrary functions. Hence the general solution of the gravitational field equations for 
collapsing strange matter in the Vaidya metric (|4|) is given by 



/ \ t-i/ s \ Mi 47ri? o q 2 (u) 

, (u, r) = F(u) + C (it) r 1/3 + — -r 3 - 

3 2r 



(20) 



Anr 2 



dF(u) + rfC(it) ^ 1/3 _ ^ dg^u) 1 



(21) 



p{u,r) 



Aitr 2 



-C(u)r 



-2/3 



q 2 (u) 
2r 2 
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(22) 
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p(u,r) 
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g 2 W 10 o 2 



(23) 



The electromagnetic current follows from the Maxwell Eq. (jlOD and is given by 
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For B = and g(w) = we obtain the solution given by Husain [11] for the null fluid pressure satisfying p = 
The energy- momentum tensor of the mixture of fluids under consideration belongs to the Type II fluids [11]. The 
energy conditions are the weak, strong and dominant energy conditions /i>0, p > 0,p > 0, p> p> and can be 
satisfied by appropriately choosing the arbitrary functions F(u) and C(u) that characterize the injection and initial 
distribution of mass and q(u) that describes the variation of the charge. The condition p > is equivalent to ^ > 
and leads to 



dF(u) dC(u) 1 /r> , .dq(u)l 
du du du r 



(25) 



imposing a simultaneous constraint on all three functions F(u), C(u) and q{u). For small values of r and for charged 
quark matter, the right hand side of Eq. ( p5| ) dominates and this energy condition could not hold. One possibility to 

satisfy Eq. ( ^5|) for all r is to assume that the function q(u) behaves so that dq d ^ — > for r — > 0. This means that 
the charge in the singular point r = is constant for all times. Alternatively, we may suppose that at extremely small 
radii matter is converted to strange quark matter so as to satisfy the energy condition. For neutral quark matter Eq. 
(^5|) is easily satisfyed by choosing > and > 0. To satisfy the condition p > for large r, we must impose 

on the function C(u) the constraint C(u) > 367r£?r 7 / 2 > 0,Vit. With this choice the condition of the non-negative 
energy density is also automatically satisfied. Due to the bag equation of state (||) we always have p > p > and thus 
the dominant energy condition holds, too. 

The radii of the apparent horizon of the metric (Q) are given by the solution of the equation 2m = r. If 
lim u _ > .oo = Fq = const., limj^oo C{u) = Co = const, and liniu^oo q(u) = qo = const, then the algebraic 

equation determining the radii of the apparent horizons is 



2F + 2C r 1 / 3 



8ttB 



(26) 



which in general may have multiple solutions. 

The singularities of the quark matter filled Vaidya space-time can be recognized from the behavior of the energy 
density and curvature scalars like e.g. R a pR a P given by 



l -C(u)r^+ q ^+A^B 
3 r 4 



(27) 



which diverges for 



III. OUTGOING RADIAL NULL GEODESICS EQUATION 

The central shell-focusing singularity (i.e. that occuring at r = 0) is naked if the radial null-geodesic equation admits 
one or more positive real root Xq [16]. In the case of the pure Vaidya space-time it has been shown that for a linear 
mass function 2m{u) — Au the singularity at r — 0,u — is naked for A < | [17]. In order to simplify calculations 
we choose some particular expressions for the functions F(u),C(u) and q(u), e.g. F(u) — au/2, C(u) = /3m 2 / 3 
and q(u) — 71*, with a > 0,(3 > and 7 > constants. With this choice the equation of the radially outgoing, 
future-directed null geodesic originating at the singularity is 

dr ! _ au _ p ^«~) 2 / 3 _ ^ ^v,y _ §7TB_ r 2 ' 
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For the geodesic tangent to be uniquely defined and to exist at the singular point r = 0, u = of Eq. (|2S]) the 
following condition must hold [16]: 
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u du , 

lim - = lim — = X . (29) 
u,i — >0 r »,i — >o ctr 

When the limit exists and Xq is real and positive, there is a future directed, non-space=like geodesic originating 
from r = 0, u = 0. In this case the singularity will be, at least, locally naked. For the null geodesic Eq. ( p8| ) condition 
(^) leads to the follosing algebraic equation: 

7A7] + a Xl + (3X1^ - X Q + 2 = 0. (30) 

With the help of the substitution X = y 3 Eq. ( p0| ) becomes a ninth order polynomial equation of the form 
f(y) = IV 9 + a V 6 + Py 5 — y 3 + 2 = 0. According to a theorem given by Poincare [18] the number of positive roots 
of a polynomial equation equals the number of changes in sign in the sequence of the non-negative coefficients of the 
polynomial g(y) = (1 + y) k f(y). The Poincare criterion indicates that Eq. ( |30| ) has two positive roots. In order to 
give an estimate of the positive roots we shall consider the expression g(y) = (1 + y) f(y) — X^io ^ k ( u ) y" ■> with 
&fc(0) = 7 > [18]. Let ^(1) denote the smallest integer for which b/. (ffc(l)) > and (ffe(l) + 1) < 0. Then we 
obtain the numbers ffc(l) and ffc(2). These numbers satisfy the relations [18] 

■ ^ <*(fc,M<£^,« = l,2, (31) 
k — Vk{s) + l k — i>k{s) 

, lim I ^ii.i = , lim ^ ( fc ' v,s)=Vs>0,s= 1, 2, (32) 

/(y s )=0,s = l,2. (33) 

Therefore it is always possible to construct a convergent sequence to obtain the positive roots of Eq. ( |30| ) . In the 
case of the neutral quark fluid, q(u) = and Eq. (|3^) is reduced to a sixth order algebraic equation also having two 
positive roots. 

IV. CONCLUSIONS 

In the present paper we considered the collapse of a strange quark fluid in Vaidya geometry. The possible occurence 
of a central naked singularity has also been investigated and it has been shown that, at least for a particular choice 
of the parameters, a naked singularity is formed. Depending on the initial distribution of density and velocity and 
on the constitutive nature of the collapsing matter either a black hole or a naked singularity is formed. The values 



of the parameters in the solution (20)-(p3[) determine which of these possibilities occurs. The solution describing the 
collapse of the quark matter is not asymptotically flat and this condition does not play any role in the formation of 
the naked singularity. Due to the presence of the bag constant B (playing the role of a cosmological constant, from a 
formal mathematical point of view) the mass function ( pp|) gives a cosmological type metric. 

Quark matter, deconfined phase of hadronic matter at high temperatures or densities may reside as a permanent 
component of neutron stars core or to form stable compact stellar objects [2,4,5]. In fact from a physical point of 
view it seems to be one of the best and more realistic candidates for the study of properties of collapsing objects. It 
also serves to ilustrate the much richer interplay that can occur among particle physics and general relativity when 
more involved quantum field theoretical models are considered. 
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